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ABSTRACT 

Hubble Space Telescope imaging of 11 fields in M31 were reduced to color-magnitude diagrams. 
The fields were chosen to sample all galactocentric radii to 50 kpc (ss 9 disk scale lengths, 
or >99.9% of the total light enclosed). The colors of the red giants at each pointing map to 
an abundance distribution with errors of order 0.1 dex in [M/H]. The abundance distributions 
are all about the same width, but show a mild gradient that flattens outside ~20 kpc. These 
distributions were weighted and summed with the aid of a new surface brightness profile fit to 
obtain an abundance distribution representative of the entirety of M31. Since we expect M31 to 
have retained most of its congenital gas and subsequently accreted material, and since the present 
day gas mass fraction is around 2%, it must be a system near chemical maturity. This "observed 
closed box" is found to suffer from a lack of metal-poor stars and metal-rich stars relative to the 
simplest closed-box model in the same way as the solar neighborhood. Models modified to include 
inhomogeneous chemical enrichment, variable yield, or infall all fit to within the uncertainties. 
As noted elsewhere, stars in the outer regions of M31 are a factor of ten more metal-rich than 
the Milky Way halo, ten times more metal-rich than the dwarf spheroidals cospatial with it, and 
more metal-rich than most of the globular clusters at the same galactocentric radius. Difficulties 
of interpretation are greatly eased if we posit that the M31 disk dominates over the halo at all 
radii out to 50 kpc. In fact, scaling from current density models of the Milky Way, one should 
not expect to see halo stars dominating over disk stars until beyond our 50 kpc limit. A corollary 
conclusion is that most published studies of the M31 "halo" are actually studies of its disk. 

Subject headings: galaxies: chemical evolution — stars: abundances — stars: giant — galaxies: 
individual: NGC 224 



1. Introduction 

We have knowledge of the spiral Andromeda galaxy (M31; NGC 224) and its elliptical companions M32 
(NGC 221) and NGC 205 that rivals and in some ways surpasses our knowledge of the Milky Way. Comparing 
the two galaxies can be a stringent test for theories of hierarchical galaxy formation. The Andromeda Galaxy 
has about the same baryonic-plus-dark mass as the Milky Way (Evans et al. 2000; Gottesman, Hunter, & 
Boonyasait 2002) but twice the number of globular clusters. Gas phase material comprises 2% of the 
baryonic total in M31, but about twice this in the Milky Way (van den Bergh 2000). Further, the Milky 
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Way is forming stars at a greater rate, with a factor of ~4 greater mass in ionised hydrogen. Most of the gas 
in M31 is concentrated in a 10-kpc "ring of fire" seen in molecular gas, atomic hydrogen, and dust, although 
OB associations appear throughout the disk, van dcr Kruit (1989) finds that the bulge of M31 contributes 
25% of the total light while Milky Way's bulge contributes more like 12%. (One should temper this result 
with the realization that the Milky Way emits more light per unit mass than does M31; in terms of mass 
the percentage contributions may be more nearly equal. Also, our own /-band bulge-disk decomposition 
discussed below yields a bulge contribution of 12%, not 25%.) 

Globular cluster systems differ in the two galaxies despite both having the same number of clusters 
per unit mass within a factor of ^2 (van den Bergh 2000). With a dividing line between metal poor and 
metal rich placed at [M/H] = —1.0, Milky Way metal rich globular clusters are concentrated toward the 
center of the galaxy. M31 has a larger number of metal-rich globular clusters than the Milky Way. As in 
the Milky Way, the metal rich M31 globular clusters lie preferentially toward the center of the galaxy, but 
unlike the Milky Way, very metal-rich globulars (of nearly solar abundance) are also found quite far from the 
center. Unlike the Milky Way, a significant subset of M31 globular clusters appear to have strongly disk-like 
kinematics (Morrison et al. 2004). 

Halo field stars in the Milky Way have an abundance distribution similar to that of the globular clusters 
spatially coresident there: with a peak at [M/H] ~ —1.5 and a broad dispersion that matches the simplest 
of chemical evolution models, the "Simple model," that assumes a closed box, full mixing, instantaneous 
recycling, and a constant heavy element yield in every generation of stars. Halo field stars in M31, if they 
do indeed belong to a halo population, are quite different. Among other authors, Harris & Harris (2001); 
Brown ct al. (2003); Rich, Mighcll, & Neill (1996): Durrell, Harris, & Pritchet (2001); Renda et al. (2005), 
and Grillmair et al. (1996) find from color-magnitude diagrams that almost all of the field stars in the outer 
regions of M31 are more metal-rich than 47 Tucanae (at [M/H]^ —0.8), the prototypical metal-rich disk 
globular cluster, with mean abundances for the giants of [M/H]^ —0.5. This is almost alarming: why should 
a galaxy of about the same mass have a halo abundance a factor of ten higher? 

The aforementioned "Simple Model" of chemical evolution is a "straw man" (easily knocked down) 
model that applies to a closed box of gas that turns into stars. However, its simplicity, along with the relative 
complexity of the alternative chemical evolution models, makes it a compelling starting point. Furthermore, 
taken as a whole, M31 inside a galactocentric radius of 50 kpc must be close to a closed box in that it does 
not seem reasonable that it permanently lost more than a pittance of its total gas mass to intergalactic space. 
Models support near-total gas retention for high-mass galaxies in non-cluster environments (Governato et 
al. 2004). Given that almost all of its gas is either still in the M31 disk (2% by mass) or has been turned 
into stars (98% by mass), then on a global level M31 is a closed box or close enough to a true closed volume 
that outflow should be negligible. 

But M31 is a closed box only in the outflow direction. Infall in M31 could have occurred throughout its 
history, perhaps episodically where each episode goes nearly to chemical completion before the next wave, 
or perhaps steadily where the M31 gas fraction has been monotonically decreasing to its minimum at the 
present epoch. In the latter case, the Simple model might be expected to apply fairly well to the M31 closed 
box without modification. In the episodic case, one would have a situation usually termed inhomogeneous 
chemical evolution. In inhomogeneous chemical evolution, different patches of the galaxy mix only slightly 
or not at all or are allowed to mix only between enrichment episodes. 

A final aspect of M31 that affects the discussion of its formation and chemical evolution is the nature of 
its stellar halo. M31 does possess a system of globular clusters that is spatially and kinematically spheroidal, 
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with a net rotation of only 80 ± 20 km s _1 (Huchra 1993). Perrett et al. (2002) finds a net rotation of 
138 ± 13 km s -1 and velocity dispersion of 156 ± 6 km s _1 for a sample of 321 clusters. Racine (1991) finds 
a globular cluster surface density profile of R~ 2 in the range 6 < i?(kpc) < 22, but a steeper drop off, or 
even a cutoff, beyond this range, van den Bergh (1969) finds that the most metal-rich globular clusters all 
have velocities of ±100 km s _1 from the projected disk velocity at their locations, indicating a more disk-like 
motion for the metal-rich end of the cluster population. More extensive data (Morrison et al. 2004) indicate 
that a significant fraction of M31 globular clusters have thin-disk kinematics, are sprinkled over the entire 
projected disk, and have abundances that span the entire metallicity range. The M31 globular clusters have 
a spread of integrated colors that covers about the same abundance range as Milky Way globular clusters, 
but contains a more equal number of clusters in each metallicity bin (Reed et al. 1994). 

The surface brightness of field stars in the outer portions of M31 as inferred from star counting follow, 
at least roughly, an exponential i? 1 / 4 profile (van den Bergh 2000). This is roughly p(R) oc i?~ 5 in the outer 
regions, as opposed to R~ 3 as one would infer for globular clusters (Racine 1991). But the continuity of 
the inner and outer parts of the minor axis (illustrated, for example, in van den Bergh (2000), figure 3.7) 
does not demand a "halo" that is separate from an "outer disk" . In fact there is no evidence of a halo at 
all from the light profile of the galaxy. Support for this viewpoint can be found from the planetary nebulae 
kinematics of Hurley-Keller et al. (2004), who say, "If M31 has a non-rotating, pressure-supported halo, we 
have yet to find it, and it must be a very minor component of the galaxy." 

In following sections we discuss HST photometric observations that sample galactocentric radii within 
50 kpc, then we present the inferred abundance distributions. (What we refer to as abundance distributions 
are called by many "MDFs," or metallicity distribution functions, and are simply the fractions of stellar 
mass observed over the range of abundance.) The observed distributions are then weighted and summed 
to simulate a globally integrated "closed box" abundance distribution for M31, which is then compared to 
solar neighborhood observations and model predictions. We conclude with a discussion of the constraints 
imposed upon the chemical evolution and formation of the Andromeda galaxy and a discussion of the relative 
dominance of disk vs. halo in the outer parts of the galaxy. 

2. Observations and Analysis 

A number of archived HST images were selected from the "planned and completed exposures" list. The 
first selection pass was simply by coordinates and whether the exposure was an image. The next step was 
to look at whether exposure times were fairly long, and if images were taken through two filters. For the 
outer fields it was clear that only the WFPC2 had a sufficiently wide field of view to count enough stars for 
meaningful results. The inner parts of M31 were observed by many programs and we sampled from those, 
ensuring that we covered as wide a range of radii as possible. Because inner fields suffer from crowding, only 
the PC chip was used there. Outer fields suffer from sparse star counts, so all available chips were utilised. 
Milky Way foreground contamination as judged from the various color-magnitude diagrams was negligible 
in all cases. 

The final fields are listed in Table 1 and illustrated in Figure 1. The table lists (1) a project-internal 
label for each field, (2) the WFPC2 filters used, (3) the semi-major axis in arcsec from our radially-dependent 
elliptical isophotal fitting, (4) the resultant galactocentric radius assuming an inclined circular disk in kpc, 
(5) the galactocentric radius in the plane of the sky (impact parameter) in kpc, (6) the original HST 
proposal number, (7) the weighting factor used for combining the different samples to make a composite 
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M31 abundance distribution, (8) the number of stars counted above the color-dependent cutoff magnitude, 
(9) the abundance median, and (fO) the abundance distribution FWHM. Many of these quantities are 
discussed further below. The HST frames were pipeline processed for basic flat fielding and photometric 
calibration. Each data set consisted of some number of frames in two filters to be coadded or "drizzled" to 
make final frames. Some of the outer fields contained additional position shifts. In these cases, each position 
was reduced separately. DAOPHOT (Stetson 1987) was used to derive fluxes for as many stars as possible in 
each frame. Zero point and aperture corrections were applied. The fluxes were corrected for charge transfer 
efficiency effects via Whitmore et al. (1999) and transformed to Johnson-Cousins V and / via Holtzman et 
al. (1995). A reddening of E(B — V) = 0.06 mag (Schlegel, Finkbeiner, & Davis 1998) was assumed for all 
fields. 

Once color-magnitude diagrams are constructed, it is a straightforward matter to overlay giant branches 
from isochrones of different abundances. Stars that lie between isochrone giant branches are assigned an 
abundance given by the mean of the flanking isochrones and binned. A faint cutoff given by M/ — 0.25(V — 
I) = —2.40 was applied so that photometric errors would not artificially broaden the inferred abundance 
distribution too much. Finally, a correction for RGB stellar lifetimes was applied to transform the observed 
star counts to mass fractions (assuming an initial mass function for low-mass stars that is constant as a 
function of heavy element abundance). More discussion about the technique can be found in Grillmair 
et al. (1996) and Worthey et al. (2004). The isochrones used were those of Worthey (1994); they have 
giant branches that fit VI cluster observations better than most existing sets over the entire range of 
abundance. Despite the faint cutoff, residual photometric error does slightly increase with width of the 
resultant abundance distribution. We ignore this effect because it is quite difficult to correctly implement 
such a deconvolution given the geometry of the color-magnitude diagram and also because rough estimates 
indicate that, with the observed widths of roughly 0.6 dex for the final abundance distributions, the increase 
of width due to photometric errors is only a few percent. 

A possible problem worth recapping is the question of age and age spreads. Age modulates giant star 
colors rather less than abundance in a logarithmic sense. As Worthey (1994) puts it, the line of null color 
change is d[log(age)]/d[log(abundance)] = —3/2 so that (for example) an age change from 8 Gyr to 10 Gyr 
can be exactly offset by an abundance change of — 2/31og(10/8) = —0.06 dex. If the underlying age spread is 
9 ±3 Gyr the corresponding abundance uncertainty is 0.08 dex. This error is slightly smaller than the one we 
infer from random photometric uncertainties (^0.1 dex). The best claim for intermediate-age subpopulations 
in the outer parts of M31 is that of Brown et al. (2003), but "intermediate age" for those authors is «6 
Gyr. A subpopulation of this age is "old" for our purposes; too old to cause any significant change in the 
abundance distribution inferred from giant colors. 

3. Abundance Results 

Figure 2 shows the resultant series of abundance distributions for the various pointings, except that the 
star-starved outer pointings have been summed. One can see a steady progression in median abundance that 
goes from metal-rich in the inner disk to relatively metal-poor in the outskirts. The widths of the distributions 
do not change much from center to edge, and are similar to the width of the abundance distribution in the 
solar neighborhood (FWHM widths are listed in Table 1, derived by averaging the top two bins to find the 
maximum, then linearly interpolating for the half- widths). 

Figure 3 shows the median abundances for each field, but split into median abundance for each chip 
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(PC, WF1, WF2, or WF3) in the cases where multiple chips were reduced. This is in order to gauge the 
effects of (poor) counting statistics. The figure shows that the outer disk-halo abundance holds relatively 
steady at [M/H] sa —0.5. Gas-phase abundances from HII regions are also plotted on the figure to show the 
present-day abundance. The nebular abundances at least roughly agree with the slope of the gradient in 
the inner 25 kpc. To the extent to which one can trust the nebular zero point, the gas is at an abundance 
roughly 0.2 dex higher than the median of the stars, which puts the gas on the high-abundance end of the 
stellar abundance distributions. This is what one would expect with Simple-model type chemical evolution 
with good mixing. With inhomogeneous chemical evolution, one would expect individual star formation 
regions to be distributed about the same way as the fossil stars. We do not emphasise this point because 
the nebular abundance zero point does have considerable uncertainty. 

The galactocentric radii of Figure 3 come from new elliptical isophotal fits of an /-band mosaic (Choi 
et al. 2002). The azimuthally-averaged surface brightness profile was extracted following the method of 
Courteau (1996) and decomposed into bulge and disk components. This is shown in Figure 4. The disk 
brightness profile is very closely exponential except for a stellar and gaseous "ring" beyond 10 kpc. The bulge 
brightness is also closer to an exponential profile than a de Vaucouleurs, with a Sersic index of n = 1.6 . 
The decomposition technique is fully described in MacArthur, Courteau, & Holtzman (2003). The surface 
brightness limit of the image was much too bright to include regions that could be considered halo regions. 
The sizes and orientations of the fitted ellipses were used to generate the Table 1 semi major axis column. 
Then, assuming an inclined circular disk and a distance of 770 kpc (Freedman & Madore 1990), these semi 
major axis values were converted to a galactocentric radius in kpc, also listed in Table 1 and used for Figure 
3. For the outer parts, we assumed an ellipticity of 0.679 and a position angle of 50°; the parameters of the 
last reliable isophote. 

The exponential isophotal fit was extrapolated to estimate the total light of the galaxy. The observed 
cumulative light profile plus the extrapolation is shown in Figure 5. Locations of HST fields are also shown 
in this figure, including the inner fields (in08, in09, and inlO) that were eventually dropped because stellar 
crowding made deriving a reliable abundance histogram impossible. The curve in Figure 5 was used to 
assign each pointing with a weight representing the fraction of the galaxy "covered" by the pointing. These 
weights are also listed in Table 1. One notices that field "In07" is heavily weighted. This is because it is the 
innermost field that was not beset with overcrowding problems. This field lies outside the bulge. This means 
that the bulge, as such, is not sampled. The fitted ratio of total bulge light to total disk light is 0.137, so 
that even if the bulge abundance distribution is very different from the rest of the galaxy it will not affect 
the all-M31 abundance distribution much. Note that the nuclear regions of M31 were studied in Worthey, 
Dorman, & Jones (1996), who conclude that M31 cannot have more than w 5% by mass contribution from 
metal-poor stars, in harmony with our results where both field In07 and the all-M31 abundance distributions 
have 3% of the mass in populations with [M/H]< — 1. 

The weights for each field were then used to construct an all-M31 abundance distribution from the 
individual sampled distributions. This global abundance distribution should be a good approximation to 
the M31 closed box. The weighted-sum abundance distribution and each of the individual distributions are 
tabulated in Table 2 (except that the outer fields are averaged). The distributions in Table 2 are normalised 
to a unit integral. The weighted sum distribution has full-width-half-maximum (FWHM) of about 0.69 dex, 
slightly broader than any of the individual abundance distributions. We compare this global abundance 
distribution to local observation and theory in the next section. 



In this notation, the exponential and de Vaucouleurs profiles have a Sersic index equal to 1 and 4, respectively. 
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4. Basic Chemical Evolution Models and Comparisons 

We can compare the global M31 abundance distribution to the solar neighborhood abundance distribu- 
tion and to simple models of chemical evolution. We do not attempt to fit the full set of radially sampled 
abundance distributions in this paper, but we recap the basics of chemical evolution and try representative 
examples of all the classical schemes (Audouze & Tinsley 1976) for stepping away from the very simplest of 
closed-box models, the Simple model. 



4.1. Analytic models 

The Simple model of chemical evolution assumes an isolated system, no infall or outflow, of one zone, 
whose total mass is constant. It begins as pure gas with a metal abundance of Z = and is well mixed at all 
times. In addition, the IMF and nuclcosynthetic yields of the primary elements in the stars remain constant. 

To recap the Simple model [c.f. Searle & Sargent (1972); Audouze & Tinsley (1976); Binney & Tremaine 
(1987)], we define Mh as the mass in gas phase made of heavy elements, M g the mass of the gas, M s the mass 
in stars and define the abundance, or, loosely, the metallicity, as Z — Mh/M g . If an incremental parcel of 
new stars is formed and only SM S remains locked in the stars, then the rest must be returned. The amount 
of gas returned by newly formed stars in heavy elements is p5M s , where p is the "nucleosynthctic yield" in 
a given generation of stars. One could assume that p is a function of time or metallicity, but the Simple 
model assumes a constant yield. We can now calculate the change in the amount of metal-rich gas, which is 
the gas returned minus the gas which formed the stars: 5Mh = p SM S — ZSM S = (p — Z)8M S . Meanwhile, 
the change in the gas metallicity is 

By combining equations and using the fact that mass is conserved such that SM S =— SM g , we obtain 

SZ = ±-[{p - Z)5M S Z8M g ] = ji-H-p + Z- Z)5M g ] = ~P^- (2) 
If the yield p is not a function of time, integration gives 

This expression for Z is a function only of the yield and the gas fraction, hence the term "Simple" for this 
scheme. Let us rework this formula to give us what we observe: the number of stars at each metallicity. If 
we let the gas fraction approach chemical completion, how many stars of what metallicity do we expect? 
The mass in stars of Z less than Z(t) at some time t is 

M s [< Z(t)\ = M g (0) - M g (t) = M fl (0)[l - exp(-^)] (4) 
so the differential fraction of mass at each Z is 



dM s _ exp(-Z/p) 
~dZ~ ~ p ' 



(5) 



- 7- 



But abundance is typically measured in terms of the logarithmic number abundance like, for example, [Fe/H], 
or [M/H] if one lumps all heavy elements together. We would like to recast equation (5) by defining F « 
[M/H] = log(Z/Z Q ). 2 Using Z = Z Q 10 F and dZ = Z dF/log(e), we arrive at 



dM = Y - , exp 
P log(e) 



P 



dF. (6) 



This distribution should be truncated at the gas fraction appropriate for the system. In the case of M31, 
the present day gas fraction is j^ 3 ^ = 0.02. 

In Figure 6 the Simple model and the global-M31 abundance distribution are compared. Note that 
M31's small present-day gas fraction means that chemical evolution in the Simple model is almost complete, 
and only a small modification in the most metal-rich bin is needed to account for the truncation of the model 
distribution. What is evident immediately is that the Simple model abundance distribution is wider, with 
an extended tail toward metal-poor stars. This is the same mismatch as the "G dwarf problem" seen in the 
solar neighborhood. 

In Figure 7 a Simple model (with a yield 0.1 dex higher than Figure 6) is compared to various solar 
neighborhood data sets. The broadest of these is Wyse & Gilmore (1995) (FWHM 0.75 dex), and the 
narrowest is J0rgensen (2000) (FWHM 0.47 dex). These widths more than bracket our derived M31 abun- 
dance distribution widths listed in Table 1 and the closed-box M31, which has a FWHM of 0.69 dex. The 
Haywood (2001) distribution is intermediate in width (FWHM 0.61 dex), but is shifted toward metal-rich 
stars. Wyse & Gilmore (1995) emphasised corrections to augment the metal-poor end of the distribution to 
account for thick disk and halo stars while Haywood (2001) emphasised selection effects against metal-rich 
stars. The differences between authors are large, illustrating that the solar neighborhood data are subject 
to a fair amount of interpretation. All of the distributions, however, are more narrow than the Simple model 
(FWHM 0.87 dex), so the solar neighborhood has a G dwarf problem for all data sets. 

Any number of modifications to the Simple model can make it fit the data better. 

The first scheme is to have a nucleosynthetic yield that starts high, then decreases with increasing abun- 
dance (an obvious way to make chemical evolution proceed more quickly in metal-poor regimes). Whether or 
not this applies in the real universe is unknown, but many authors have been driven to consider the idea of 
a variable initial mass function that might drive a variable yield (Bresolin, Kcnnicutt, & Garnett 1999; Cen 
2003a, b; Scannapieco, Schneider, & Ferrara 2003; Schneider et al. 2002; Chabrier 2003; Larson 1998; Padoan 
et al. 1997; Chiappini et al. 2000; Chiosi et al. 1998; Kroupa & Weidner 2003). Alternatively, supernovae 
themselves may have a different character at lower abundance. 

We modify the Simple model to have a yield that begins high and then decreases with increasing Z 
as p — tf^- We have not seen this particular variant in the literature; we call it the "rational decreasing 
yield model" , but really it is not a physical model and serves only to illustrate how easily one can fit 
the observations. This scheme is a three-parameter model with the gas fraction remaining as the primary 
parameter, but the yield now parameterised with po and e. Substituted into equation (2) and integrated, 



2 Note that this is not precisely correct. If "heavy" elements are those other than hydrogen or helium, then one should 

Z I (\ z) 

account for the fact that very metal-rich stars deplete their light elements and define F = log % . /\i — z .) ' H° wever i if one 
considers stars from very metal-poor to a few times the solar abundance of Zq Kb 0.02, one can safely retain the approximate 
form. 



- 8 - 



the expression for p yields a quadratic with positive root 



Z(t) = -e+ (e 2 - 2p ln 
The analog of equation (4) for this model is 



M g (t) 



M 9 (0) 



) 1/2 . 



M s [< Z(t)} = M g {0)[ 1-exp 



(Z + e) 2 -e 2 



The analog of equation (5) is 



dM s Z + e 



exp 



dZ po 

And converting to ten-based logarithmic units yields 



2po 

(Z + ef + e 2 - 



2 Po 



JJ1# Z Q 10 F + e 

dM s = — exp 

Po 



(Z e 10 F + e) 2 - e 2 



2 Po 



Z n 10 F 



log(e) 



dF. 



(7) 



(8) 



(9) 



(10) 



Figure 8 illustrates the improvement in fit that can be achieved with a variable yield. The parameter p 
controls the abundance of the peak of the distribution, while e controls the width. The scale on the y axis 
and the thin decreasing line represent the variable yield, computed with p — 0.00019 and e = 0.004. The 
shape of the resultant abundance distribution is not sensitive to modest changes in the parameters; a good 
feature of this modified model. Another good feature is the plateau at the very metal-poor end, but there 
are mild mismatches at various points as well. The width is FWHM = 0.62 dex for the parameters listed. 



4.2. Inhomogeneous Enrichment Models 

Suppose that, in a cube of gas, enrichment occurs in patches scattered throughout the cube. These 
patches may or may not overlap spatially, so that if little mixing occurs then a variety of abundances will 
coexist in the gas. Such inhomogeneous chemical evolution schemes been studied in several works [e.g. 
Tinsley (1975), Searle (1977), and Malinie et al. (1993) ]. 

First, we compare our M31 composite abundance distribution to a model with zero mixing (Oey 2000). 
This model has five parameters: the volume filling fraction Q of the enriched bubbles, the number n of 
enrichment generations, the mean a and standard deviation a of the abundance distribution typical of a 
collection of local star formation events, and the present day gas fraction /x = M g (t)/M g (0). We assume 
that each generation of star formation is described by a Gaussian abundance distribution; i.e., we stop at 
Oey's equation (8) rather than going on through equation (13). 

Some Oey (2000) models are plotted in Figure 9, with parameters listed in the figure itself. The broad 
model (FWHM 0.76 dex) with n = 400 and Q = 0.8 is similar to what is plotted in her Figure 4 that is 
meant to match the solar neighborhood data, although Oey uses a different gas fraction. The other, narrower 
distribution (FWHM 0.55 dex) has the number of events multiplied by the volume filling factor of order one 
(nQ ~ 1). We worry that this corner of parameter space is not ideally suited for spiral galaxies with long 
star formation histories, but the model does fit quite nicely. 

We also compare the M31 abundance distribution to the inhomogeneous chemical evolution model of 
Malinie et al. (1993) in which mixing is allowed between, but not during, inhomogeneous star formation 
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events. This model was built to emulate the solar neighborhood situation and, like the Oey model, has 
5 parameters; the number of mixing events N, the fraction of gas consumed in each event /, an initial 
metallicity Zq, a metallicity dispersion SZ, and a yield y. Figure 8 illustrates the match of the Malinic ct al. 
(1993) model (dashed line) if the initial metallicity Z = and the fraction of gas consumed in each event 
(/) is set so that the final gas fraction is 0.02 instead of 0.2 as in their paper. Other parameters are left as in 
Malinie et al. (1993) (N = 100, 5Z = 0.01, and y = 0.72Z Q ). Setting the number of events N equal to 100 
implies that there was quite a lot of mixing in the history of M31. The Malinic model (FWHM 0.63 dex) is 
quite successful at roughly matching the M31 closed-box abundance distribution (FWHM 0.69). The only 
parameter changes from the Milky Way case were the expected ones: an adjustment for a lower gas fraction 
in M31 and an initial metallicity of zero that is appropriate for a closed box situation. 



4.3. Infall 

We can also massage the shape of the abundance distribution by feeding gas into the box from outside, 
even as chemical evolution proceeds within the box. This could be done with complete mixing or with 
inhomogeneity. We illustrate the former option using the parameterization of Harris & Harris (2002). This is 
a five-parameter model: the metallicity of the infalling gas, a constant yield, and a constant gas consumption 
rate, expressed as the fraction of gas consumed per time step, e. The rate of infall is constant for some time 
Ti, then decays exponentially as exp(— i/r 2 ), where t is an integral number of time steps. 

Figure 10 shows representative models compared to the M31 data. The most sensitive input parameters 
are the yield and the Ti infall timescale. In Figure 10 a model that fits the data (with FWHM=0.68) is 
plotted along with models in which the decay timescale is varied by a factor of 1/4 and 4. A more constant 
star formation rate approaches a delta function shape in the abundance distribution, while a very short star 
formation time makes the distribution very broad. While the yield might be a universal constant, infall 
timescales are not, so this model gives a prediction that galaxies with short formation times or that formed 
in bursts (like at least some elliptical galaxies) should have much broader abundance distributions than 
galaxies with long formation timescales (Sc or Sd or LSB galaxies) unless infall rate is directly proportional 
to star formation rate; the model abundance distribution is invariant in this case. Due to the sensitivity of 
the model to the infall timescale, the infall rate would have to be proportional to the star formation rate to 
a high degree of accuracy across Hubble types and across quiescent or bursty star formation modes in order 
for it to apply to real galaxies, which tend to have very similar abundance distributions: Worthey, Dorman, 
& Jones (1996) note that, so far, every galaxy studied that is bigger than M32 seems to suffer a "G dwarf 
problem" of having a narrow abundance distribution compared to the Simple model. [Objects include the 
Milky Way, M32, M31, and the nuclear regions of elliptical galaxies. The outer regions of disturbed giant 
elliptical NGC 5128 have also been shown to have a narrow abundance distribution with a FWHM ~ 0.65 
dex (Harris & Harris 2002; Rejkuba et al. 2004, 2005)]. This argues for a universal phenomenon rather than 
a special case scenario. Observed star formation rates in gas-rich galaxies range from < 10~ 2 M Q yr _1 in 
Sc/Sd/LSB spirals (Burkholder et al. 2001) to > 10 3 M yr 1 in ultraluminous infrared galaxies (Smail et 
al. 2002). One would have to require that infall rate and star formation rate are synchronized to about a 
factor of two over this range in order for this exact infall model to apply. On the other hand, infall must 
be an important component of galactic evolution. Plausibly, alternative methods of parameterization may 
show less volatility. 
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5. 



Discussion 



5.1. 



Closed Box Models 



In summary, the comparisons of the observed M31 closed box abundance distribution with chemical 
evolution models showed that, while the Simple model is too broad, all kinds of slightly modified chemical 
evolution schemes fit the narrower observed abundance distribution of closed-box M31 stars. Combined 
with the result of Worthey, Dorman, & Jones (1996) that elliptical (i.e. spheroid-dominated) galaxies also 
have a narrow abundance distribution compared to the Simple model, the "narrowness" of the abundance 
distribution appears to be a universal phenomenon in all large galaxies, not just for disk galaxies. This 
is circumstantial evidence that the abundance distribution is generically narrower than the Simple model 
suggests. 

Since outflow is ruled out for closed boxes, the remaining classical solutions to this problem (Audouze & 
Tinsley 1976) are (1) prompt initial enrichment, (2) infall without outflow, (3) dropping the assumption of 
instantaneous recycling, i.e. allowing for time-delayed yields and allowing inhomogeneous enrichment. Our 
Z-depcndent yield is an implementation of (1), and it works fairly well except that the parameterisation 
is not physical. Scheme (2), infall, easily matches almost any abundance distribution with suitable fine- 
tuning of yield, mixing, and infall parameters, as we show with the full-mixing Harris & Harris (2002) 
parameterization. The problem with infall models as applied to galactic closed boxes is that the output 
abundance distribution is very sensitive to the infall timescale, so that the infall timescale would have to be 
nearly the same for ellipticals as it is for spirals since they have very similar abundance distributions. This 
objection is nullified if the infall timescale is always directly proportional to the star formation rate to about 
a factor of two in all star formation environments. Scheme (3) was partially explored with inhomogeneous 
models. The easiest match included substantial mixing of gas during chemical evolution. Delayed injection 
of heavy elements serves to add metal-poor stars to the system, and so would make the fit worse if it were 
included. We are therefore left with a theme often voiced about chemical evolution; that there is a uniqueness 
problem. Several models with different assumptions can match the data quite well (Tosi 1988). 



As a prelude to further discussion, we ask "at what radius do we expect M31 halo stars to dominate 
over disk stars?" assuming that the M31 halo resembles that of the Milky Way. Disks are usually modelled 
with an exponential profile of the form I — I exp(— R/h), where h is the scale length and / is a surface 
brightness, not a density. The Milky Way halo is found to follow a power law density p — poR~ 3 5 outside 
of 5 kpc (van den Bergh 2000). Provided we can trust these functional forms, an exponential will fade faster 
than a power law, so a halo will always dominate outside some radius. 

To evaluate the surface density of the halo, normalised to the local surface density of the disk, we use 
the results of Sandage (1987), integrating his Table 1 (kinematically derived densities as a function of height 
above the disk) and adopting local density ratios for disk:thick disk:halo of 500:30:1 (and for our purposes 
we simply subsume the thick disk into the disk). The problem is then fully specified except for the disk 
scale length h. Figure 11 shows the prediction for the F-band Milky Way scale length h = 4.2 kpc (van den 
Bergh 2000) . Van den Bergh notes that estimates of 2 to 6 kpc for the Milky Way scale length h exist in the 
literature, part of this spread being a function of the wavelength studied. The Chen et al. (2001) models are 
also shown in Figure 11 exactly as integrated from the published density model (including a shallower halo 



5.2. 



Disk Dominance 
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density power law of —2.5 and disk scale length of h = 2.25 kpc). These results indicate that, in the Milky 
Way, halo stars become more common than disk stars at galactocentric radii of 35 kpc < R < 45 kpc. 

The scale length of M31 is more secure than that of the Milky Way. van den Bcrgh (2000) quotes h = 5.7 
kpc for the F-band, and our own /-band decomposition gives h = 5.6 kpc [or 5.7 kpc if the slightly longer 
distance of 784 kpc (Stanek & Garnavich 1998) is used]. For a scale length of this size, and also doubling the 
density of the M31 halo because M31 has twice the number of globular clusters, Figure 11 indicates that we 
should not expect to see the halo dominate at any radius that we sample, and that halo stars should become 
relatively numerous only outside a radius of 50 kpc from the center of M31. 

This should be slightly tempered by a geometrical consideration: If the halo is spherical, a line of sight 
intersects the halo at an impact parameter Rh ~ DO, where D is distance and 9 is the angle between the 
line of sight and the M31 nucleus. The impact parameter is the same as the galactocentric radius in the 
inclined disk (Rd) only along the major axis of the isophote. For most lines of sight, the projected disk is 
foreshortened so that Rh < Rd, and the foreshortening is at a maximum along the minor axis of an isophote 
(Walterbos & Kennicutt 1987). For gauging this possible amplification effect, we list both radii in Table 
1 [the amplification of the halo is roughly (Rd/ Rh) 2 ' 5 assuming the halo density profile is p cx R 3 - 5 ]. We 
find that fields near the minor axis can have the halo boosted by a factor of order ten from this effect. Of 
course, if the halo is flattened instead of spherical (as the Milky Way's is believed to be) the boost rapidly 
diminishes. And there is evidence from globular clusters that the halo density falls off faster than p oc R~ 3 5 
outside of 30 kpc (Racine 1991). It is possible that, with the exception of nearly edge-on spirals where 
projection effects dominate, the disk stars outnumber the halo stars at all radii. We note that the metallicity 
remains about [M/H] ~ —0.5 even in fields as far as Rh = 30 kpc (Durrell, Harris, & Pritchet 2004, 2001). 

If we are seeing disk stars in all fields, then explaining the relatively metal- rich mean abundance of 
[M/HJrs —0.5 becomes quite natural compared to trying to explain why the halo of M31 is a factor of ten 
more metal-rich than the Milky Way halo. If we are seeing disk stars, then the M31 halo can be presumed 
to be alive and well and very similar to the Milky Way halo, which, after all, is supposed to compose less 
than 2% of the mass of the Galaxy. Alive and well, but, as yet, unseen except for the metal-poor-selected 
stars of Rcitzel & Guhathakurta (2002). 

Many authors have called the outer regions of M31 "halo," but for the most part this is a term of 
convenience, and not indicative of true pressure-supported halo status since we are almost completely ignorant 
of the kinematics of these stars as of this writing. Parantheticallly, our ignorance is about to be relieved 
as large telescopes are beginning to measure many velocities of giant stars in the uncrowded outer regions 
of M31 [cf. Reitzcl & Guhathakurta (2002)]. It is interesting to wonder how interpretations change if 
disk status is assumed. Brown et al. (2003) find an "intermediate age" subpopulation in a field at 11 kpc 
galactocentric radius in the plane of the sky (we derive 11.4 kpc from the coordinates given in their paper 
plus our assumed distance), or 35.1 kpc galactocentric radius on the extrapolated M31 disk. If this is disk, 
then finding non-ancient stars seems perfectly normal, and we may wonder why there are no younger stars 
present. Perhaps this is evidence that star formation in the outer M31 disk was quenched after about half a 
Hubble time. Preliminary Keck/DEIMOS velocities of stars around the the Brown et al. (2003) field indicate 
a velocity dispersion of 85 km s _1 centered on the systemic velocity of M31 (J. S. Kalirai (2005), private 
communication). Compared to the 156 km s _1 line-of-sight velocity dispersion of the M31 globular clusters, 
this is neither a dynamically cold thin disk nor halo, but probably something akin to the Milky Way's thick 
disk component and very much in line with our speculation of a dynamically heated ancestral disk. Recently, 
Chapman et al. (2005) report disklike kinematics between 15 and 40 kpc, but with a much smaller velocity 
dispersion (30 km s^ 1 ). 



- 12 - 



By counting red giants in the vicinity of M3 1 , Ferguson et al. (2002) established the presence of dynamical 
clumping or streaming at radii out to 55 kpc in an elliptical survey region shaped like an extrapolation of 
the disk. Some substructure detected strongly resembles merging events, but some of the features inside 
^30 kpc also look like flocculent spiral structure. Companions M32 and NGC 205 project on the sky inside 
these structures. The portions of the Ferguson ct al. (2002) Fig. 2 map outside ^30 kpc look rather uniform 
(except for the substructures) but this may be due to Galactic foreground stars that were not subtracted. 
The Ferguson et al. (2002) maps certainly do not rule out that the outer parts of M31 are dominated by a 
disturbed and tidally shredded disk. 

Our synthesis of these considerations is that M31 formed a gaseous, well-mixed disk and a Milky Way- 
like halo in the first half of a Hubble time. This ancestral disk need not be oriented exactly the same way 
as the present-day disk. Chemical enrichment proceeded in a "disky" way via rotational mixing, galactic 
fountains, and the occasional supernova wind. Over time, the gas supply ran short so the outer disk stopped 
forming stars. Additionally, late accreting satellite galaxies made their dynamical presence felt, adding stellar 
streams and tidally disturbing the disk stars that were already present. The advent of more measurements 
of stellar kinematics in the outer disk of M31 should reveal the presence of streams from disrupted satellites 
and also the structure of the postulated disturbed ancestral disk that may have a variable velocity dispersion 
as a function of galactocentric radius or even gaps and warps caused by strong perturbers. 

In particular, if M32 has been orbiting M31 for more than a few orbits it will certainly be well on the 
way to clearing a gap in the ancestral stellar disk due to tidal resonance effects. Finding the size and shape 
of such a gap could be used to constrain M32's orbit and infer how long it has been a satellite of M31. 

We thank P. Guhathakurta and P. I. Choi for providing M31 image mosaic data. P. B. Stetson provided 
both the DAOPHOT program and valuable advice. S. Oey provided feedback of inestimable worth. S. C. 
and L. A. M. acknowledge support from the National Science and Engineering Research Council of Canada. 
We thank J. S. Kalirai for permission to publish his preliminary value for the velocity dispersion at Rd = 35 
kpc, and J. L. Serven for coining the term "Worthey Gap," although we are not sure if it applies to M31 or 
G. W. This work was supported by grant AR-08745.01-A from Space Telescope Science Institute. 
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Fig. 1. — Program field locations, labeled by our internal ID, are superimposed on M31 isophotes by Hodge 
& Kennicutt (1982). The three "out" fields happen to be globular cluster fields (with the cluster subtracted) 
and the cluster names are in parentheses. Isophotes for /j,b = 25, 24, and near-nuclear 21 mag arcsec -2 arc 
shown. A rough guess for [ib = 27 is sketched. NGC 205 is marked by a dot located inside the /i = 25 
"finger" north of the M31 nucleus, and M32 is marked by a dot just to the left of the "mil" label. 
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Fig. 2. — Derived abundance distributions for each field, except that the outer fields have been averaged. 
Each distribution has been normalised to unit integral, bins are 0.1 dex wide, and each vertical tick mark is 
0.1. Labels for each field appear to the right, along with the isophotal ellipse semi-major axis value associated 
with each field. A mild abundance gradient is clearly seen. 
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Fig. 3. — Median abundance for each field plotted as a function of galactoccntric radius. Multiple points 
appearing at one radius indicate separate results for individual WFPC2 chips. Results from HII regions from 
Blair et al. (1982) (squares) and Dennefeld & Kunth (1981) (triangles) are also shown. 
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Fig. 4. — Azimuthally-averaged surface brightness profile from elliptical isophotal fits to the Choi et al. (2002) 
/-band image. The bulge light was fit with a Sersic model of index n — 1.6, effective radius, r e = 234'.'4 
(=0.87 kpc), and effective J-band surface brightness /j, e — 16.6 mag arcsec -2 . The disk is described with an 
exponential profile of scale length h = 1501'.'6 (=5.6 kpc) and a central I— band surface brightness /j,q = 17.6 
mag arcsec -2 . The dotted lines show the surface brightness profile variations due to a 1% sky error. Residuals 
from the fit are shown in the lower panel. The x 2 estimators and decomposition technique are fully described 
in MacArthur, Courteau, & Holtzman (2003). 
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Fig. 5. — The measured cumulative surface brightness profile (inside 6400"; solid line ) given as a fraction 
of the total extrapolated flux and the assumed extrapolation (outside 6400"; dotted line) are shown as a 
function of the semi-major axes of fitted elliptical isophotes. Locations of the target fields are also shown. 
Under the assumption of a circular disk geometry, galactocentric radii are given on the top scale assuming a 
distance of 770 kpc for M31. The three innermost fields were dropped from the analysis due to severe stellar 
crowding that prevented accurate photometry. 
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Fig. 6. — The global, weighted-sum M31 abundance distribution is compared to the Simple model with yield 
[M/H] p «log[p/Z Q ] = -0.2 and a gas fraction of 2%. 
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Fig. 7. — A Simple model with yield [M/H] p = —0.1 and gas fraction of 2% is compared to three independent 
literature abundance distributions derived for the solar cylinder. The literature references are Wyse & 
Gilmore (1995), J0rgensen (2000), and Haywood (2001). 
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Fig. 8. — The M31 abundance distribution is compared to the rational decreasing yield model with p = 
0.00019 and e = 0.004 (dotted line). The thin line and the y-axis scale illustrate the yield p computed with 
these assumptions. The slightly ragged dashed line is a Malinie et al. (1993) model with parameters as in 
their paper except that the initial metallicity is zero and the / parameter has been set to 0.075 so that the 
final gas fraction is 0.02. 
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Fig. 9. — The M31 abundance distribution is compared to two inhomogeneous models computed using the 
formalism of Ocy (2000) with \i = 0.02 and other parameters as noted in the plot. 
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Fig. 10. — The M31 abundance distribution (thin line) compared to infall models computed using the 
parameterization of Harris & Harris (2002) with yield = 0.7Z Q , fraction of gas consumed per time step 
e = 0.02, and n = 10 time steps of constant infall before exponential decline sets in. The subsequent decay 
time constant is t 2 = 65 time steps for the model that matches the observation (bold), and r 2 is set to four 
times longer for the narrow model whose peak is off the graph (at 0.14; dashed), and four times shorter for 
the broad model (dotted). 
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Fig. 11. — The ratio of halo to disk surface density (that is, projected on the x, y plane) for some models of 
the Milky Way. Values greater than one indicate that most stars seen by an observer distant from us and 
looking down on the disk would be halo stars, and values less than one indicate that disk stars dominate the 
light. The Chen et al. (2001) model has a scale length of 2.25 kpc, and a model based on Sandage (1987) 
plus van den Bergh (2000) 's preferred scale length of 4.2 kpc is also shown. The Milky Way in these models 
has a halo that becomes visible somewhere between 35 and 45 kpc. Both of these models arc modified to 
emulate M31 by increasing the disk scale length to 5.7 kpc and (conservatively) doubling the density of the 
halo. The figure shows that, with a Milky Way-style halo, we should expect to see very few halo stars at 50 
kpc from the center of M31. 



